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YURI A. KORDYUKOV AND ANDREY A. YAKOVLEV 

Abstract. We prove an asymptotic formula for the mimber of integer 
points in a family of bounded domains with smooth boundary in the 
Euclidean space, which remain unchanged along some linear subspace 
and expand in the directions, orthogonal to this subspace. A more 
precise estimate for the remainder is obtained in the case when the 
domains are strictly convex. 

Using these results, we improved the remainder estimate in the adi- 
abatic limit formula (due to the first author) for the eigenvalue distri- 
bution function of the Laplace operator associated with a bundle-like 
metric on a compact manifold equipped with a Riemannian foliation in 
the particular case when the foliation is a linear foliation on the torus 
and the metric is the standard Euclidean metric on the torus. 



1. Statement of the problem and the main results 

A classical problem on integer points distribution consists in the study of 
the asymptotic behavior of the number of points of the integer lattice in a 
family of homothetic domains in . This problem is originated in the Gauss 
problem on the number of integer points in the disk, where it is directly 
related with the arithmetic problem on the number of representations of 
an integer as a sum of two squares, and sufficiently well studied (see, for 
instance, books O El [HIS] and the references therein). 

In this paper we investigate much less studied problem on counting of 
integer points in a family of anisotropically expanding domains. More pre- 
cisely, let F be a p-dimensional linear subspace of M" and H = F-^ the 
g-dimensional orthogonal complement of F with respect to the standard 
inner product (•, •) in M", p + q = n. For any e > 0, consider the linear 
transformation : M" —?■ M" given as follows: 

jx, if X G F, 
\ e X, if X G H. 
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For any bounded domain S with smooth boundary in M", we put 

(1.1) n,(5) = #(r,(5)nz-), e>0. 

The main goal of the paper is to study the asymptotic behavior of the 
function 12^(3) as e ^ 0. Before we state the main results of the paper, we 
introduce some auxiliary notions. 

Let r = n -F. r is a free abelian group. Denote by r = rankF < p the 
rank of F. For r > 1, denote by {ii,i2, ■■■tr) some base in F. Let V be the 
r-dimensional subspace of spanned by (£i,£2; ■■■tr). Observe that F is a 
lattice in V . Denote by Q C ^ the parallelepiped spanned by the vectors 
{ii,i2, ■■■ir) and by IQI its r-dimensional Euclidean volume: 

|g| =vol,(4,^2,-4) =vol(y/F). 

Let F* denote the lattice in V, dual to the lattice F: 

F* = {7* G y : (7*,F) C Z}. 

For r = 0, the groups F and F* are trivial, and it is natural to put \Q\ = 1. 

For any x G V, we denote by Px the (n — r)-dimensional afhne subspace 
of R", passing through x orthogonal to V. 

Theorem 1.1. For any bounded open set S with smooth boundary in M", 
the formula holds: 

(1.2) n,(5) = ^ ^ voln-r-(Pr nS) + 0(ei^^''), e ^ 0. 

Note that, in the general case, the leading term in the asymptotic formula 
for ne(5') as e — >■ was unknown. 

Remark 1. In the case when F is trivial, we have p = r = 0, q = n. The 
problem is reduced to the classical problem on the asymptotics of the number 
of integer points in a family of homothetic domains in M"'. The formula (jl.2jl 
is reduced to the classical formula, going back to Gauss: 

n^(5) = e-"vol„(5) + 0(ei""), e ^ 0. 

Remark 2. Li the case when F is trivial, we have r = 0. The formula (|1.2p 
takes the form 

In order to obtain a more precise estimate for the remainder, we need to 
impose some restrictions to the boundary of the domain S. We give just 
one result of similar type. 

Theorem 1.2. For any bounded open set S with smooth boundary in M" 
such that, for any x £ F, the intersection Sri{x + H} is strictly convex, the 
formula holds: 

(1.3) n,(5) = ^ J] vol„_,(P^. n5) + 0(e'^-«), e ^ 0, 
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where 



2(p-r+i) if^<P-r 
^ tf^>p-r. 



L g+l+2(p-r) 

Remark 3. In the case when F is trivial, we get k = 2n/{n + 1), and the 
formula ()1.2p is reduced to the following formula: 

n,{S) = e-"vol„(S) + 0(e~"+2"^), e ^ 0. 
This formula was proved by Randol [12l [13] . 

Example 1.3. As the simplest non-trivial example, we consider the case 
when n = 2 and p = 1. Thus, let F be the one-dimensional linear subspace 
of spanned by (1, a) € M?. Then its orthogonal complement H is spanned 
by {—a, 1) € M^. One can distinguish two cases. 

Case 1: a Q. In this case, F is trivial, therefore, r = 0. Moreover, V* 
is trivial, and, as mentioned above, it is natural to put \Q\ = 1. For any 
bounded domain S with smooth boundary in M^, we get 

n,[S) = e-^area(5) + 0(e"^/^), e ^ 0. 

Case 2: a G Q. Write a = |, where p and q are coprime numbers. Then 
F is generated by = {q,p). Therefore, 



\Q\ = \h\ = vV^+^- 

The dual lattice F* is generated by (i- An arbitrary element of F* has 
the form 

7* = „ ^ M l, keZ. 
pz _j_ qz 

The corresponding subspace P^* is the line L^. on the plane M? given by 
the equation qx + py — k = 0. Therefore, for any bounded domain S with 
smooth boundary in M^, we get 

n,(5) = e'i^=L=V|5nL,| + 0(l), 8^0. 

The problem on counting of integer points in anisotropically expanding 
domains in a slightly different context was studied in considerable detail in 
[m El [ini [H] . More precisely, these papers were devoted to the estimates 
of the number A^(S', F) of points of a lattice F in M" lying in a bounded 
domain S for special domains and lattices. Let us briefly describe some 
results of these papers and show how they can be applied to the problem 
under consideration. 

Let A: be a totally real algebraic number field of degree n, a the canonical 
embedding of k in the Euclidean space M", M C k an arbitrary Z-module 
of rank n and F^ = <7(M) the corresponding algebraic lattice in M". Let 
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n C M"' be an n-dimensional parallelepiped centered at the origin and with 
edges parallel to the coordinated axes: 

n 

n = n(-«i'«i)- 

For A = (Ai, ...,Xn) € M."", we put 

n 

NmA = JJ Aj. 

i=i 

Let {ei, . . . , Sn} be the standard basis of M". Any vector A = (Ai, . . . , An) G 
M" defines a linear transformation of M" by 

A • = AjCj, jG{l,...,n}. 

By [m Theorem 1.1] (see also [Ml [TTj), for any A G M" we have the 
estimate 

<C[ln(2 + |NmA|)]"-\ 

where d(TM) is the volume of a fundamental domain of the lattice rj\,/ and 
C > is a constant, independent of A. 

In [ini IH] , this result was extended to the case of algebraic lattices asso- 
ciated with an arbitrary algebraic number field as well as to a wider class of 
domains. As an application, the authors obtain nontrivial remainder esti- 
mates in the asymptotic formula due to Lang [9] for the number of elements 
of an algebraic number field in a parallelotope determined by the canonical 
system of valuations. 

The study carried out in [Ml [I5l [lOl [H] shows that the remainder estimate 
in the asymptotic formula for N{X ■ S, T) turns out to be very sensitive to 
the Diophantian properties of the lattice T and to the geometric properties 
of the domain S. 

In order to apply the results described above to the problems considered 
in this paper, let us fix p G {1, . . . , n} and make use of the formula (jl.4p for 
A = A(e) G M" of the form 

^ [e-^, ifp+l<j<n, 

for some e > 0. We get 

<C7i|lne["-\ e^O, 

where Ci > is a constant, independent of e, and p + q = n. 

On the other hand, the lattice Tm can be represented as Tm = A['L^) 
with some linear isomorphism A of M". Therefore, we have 

iv(n,rM) = A^(n,^(z")) = N{A~\ii),ii'). 



(1.4) 



iv(A-n,r 



NmA| 



voin(n) 



d{TM) 



N{\{e)-n,TM)-e- 



, V0ln(n) 

d{TM) 
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It is not difficult to see from fiere tliat 

N{\{e)-U,rM) = n,{A~\^)), 

wliere the riglit-hand side of tliis identity is defined by (jl.ip with the sub- 
space F C M" spanned by ^"-"^(ei), . . . , A~^{ep) and the subspace C M" 
spanned by ^-^(ep+i), . . . ,^-^(e„). Observe that Z"" n F = {0}. 
Thus, finally we get 

n,{A-\U)) = e-'^voUA''m + 0(1 Inel"-'), e ^ 0, 

which is the formula (|1.2|) with a more precise remainder estimate in the 
case when S = A~^U is a parallelepiped centered at the origin and with 
edges parallel to the vectors A~^{ei), . . . ,^~^(e„). Note that, in this case, 
the subspaces F and H have rather special form and, in general, are not 
orthogonal. 

Similarly, one can use the results of [lOtlllj and get more precise remainder 
estimates in (jl.2p for some special subspaces F and H and domains S. 

It would be interesting to continue the study of remainder estimates in 
the asymptotic formula (ll.2p . depending on geometry of a domain S and 
properties of F and H. 

The authors are grateful to the referee for useful remarks. 

2. Applications to adiabatic limits 

It is well known that the Gauss problem on counting of integer points 
in the disk is equivalent to the problem on the asymptotic behavior of the 
eigenvalue distribution function of some elliptic differential operator on a 
compact manifold, namely, of the Laplace operator on a torus. In the case 
under consideration, there is also an equivalent asymptotic spectral prob- 
lem, namely, the problem on the asymptotic behavior of the eigenvalue dis- 
tribution function of the Laplace operator on a torus in the adiabatic limit 
associated with a linear foliation. 

Consider the n-dimensional torus T" = M^/Z". Let be a linear foliation 
on T*^: the leaf of T through x S T" has the form: 

L^ = x + F mod Z". 

The decomposition of M"" into the direct sum of subspaces = F ^ H 
induces the corresponding decomposition g = Qf + Qh of the Euclidean 
metric into the sum of the tangential and transversal components. Define a 
one-parameter family of Euclidean metrics on M" by 

ge = gp + e~'^gH, e>o. 

We will also consider the metrics gs as Riemannian metrics on T". 

Let A = (ai, . . . , o„) G M" be some point. Define a 1-form A on T" by 

n 

A = ajdxj. 
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Consider the operator d — 27riA, acting from C°°(T") to the space Q^{T'"') 
of smooth 1-forms on T", where d is the de Rham differential, and A is the 
multiphcation operator by A. Let {d - 27riA)*^ : 0^(T") C°°(T") be 
the adjoint of d — liriA with respect to the inner products in C°°(T"') and 
0^(T") determined by gs- 

For any e > 0, consider the operator in C°°(T") defined by 

He = {d- 2iTiA)*g^{d - 27riA). 

In the local coordinates (xi, X2, • • • , Xn) of the space M", the operator is 
written in the form 




where gi are the elements of the inverse matrix of g^ ■ 

The operator H,. can be considered as the magnetic Schrodinger operator 
on the torus T", associated with the metric g^, with the constant magnetic 
potential A. It has a complete orthogonal systems of eigenfunctions 

Uk{x) = e'^'''^'''''\ xGM'^, A:GZ", 

with the corresponding eigenvalues 

n 

Xk = (2^)2 life - AW^., = (27r)2 ^ 5f {kj - aj) {h - a,) . 

Denote by A'^j(A) the eigenvalue distribution function of H^: 
Ne{X) = ^{k e : Afc < A}, AG M. 
It is easy to see that 

n,{B^{A)) = N,{4Tr^X), A G M. 

Thus, the problem on the asymptotic behavior of the number n£{B^(A)) 
of integer points in the ellipsoid T^{B^{A)) as e — )• is equivalent to the 
problem on the asymptotic behavior of the eigenvalue distribution function 
Ni,{X) as e — 7> 0. The limiting procedure e — > will be called adiabatic 
limit. This notion was introduced by Witten in 1985 in the study of global 
anomalies in string theory. We refer the reader to a survey paper [7J for 
some historic remarks and references. 

In [5] (see also [6]), the first author computed the leading term of the 
asymptotics of the eigenvalue distribution function of the Laplace operator 
associated with a bundle-like metric on a compact manifold equipped with 
a Riemannian foliation, in adiabatic limit. The linear foliation on the torus 
is a Riemannian foliation, and a Euclidean metric on the torus is bundle- 
like. As a straightforward consequence of Theorem 11.21 obtain a more 
precise estimate for the remainder in the asymptotic formula of ^ for this 
particular case (see also [16]). 
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Theorem 2.1. For A > 0, the following asymptotic formula holds as e ^ 0: 

,, / \ \{n-r)/2 

where uJn-r is the volume of the unit ball in R"^*" and 

[ g+l+2(p-r)' ^f—>P-f- 

3. Proofs of the main results 

We will use a method based on the Poisson summation formula and the 
method of stationary phase (Van der Corput, Randol [121 [13]) Colin de 
Verdiere [1]). First of all, we observe that we have the inclusion 

(3.1) z"c y P^*. 

Indeed, let k ^ 17^. Denote by ny : R" V the orthogonal projection on 
V. Then, for any 7 G F, we have 

{7Tv{k),-/) = (A:,7) GZ, 

since 7 G and j (z F C V . Hence, Try (A;) G F*, that immediately implies 

For any 7* G F*, denote 

Z!;, = Z"p|P^. = {A: G Z" : nvik) = 7*}. 

We identify the afhne subspace P^* with the linear space V-^, fixing an 
arbitrary point k-y* G Z"* : 



It is easy to see that 
where 



l"], = k^, + F^ 



F-L = Z" 

is a lattice in V-^. 

Let us observe the following relation, which will be needed later: 

(3.2) vol(y^/F^) = \Q\. 

For its proof, let us choose some base (^1, ...,£,.) in F. Denote by (£|, . . . , £*) 
the dual base in F*: {£i,£*) = 5ij for any i,j = 1, . . . ,r. For any i = 1, . . . ,r, 
choose some k* G Z" such that 7ry(A;*) = £*. Let [k^, . . . , k;^_^) be a base in 
P-*-. Using (I3.1j) . it easy to show that the system (fc*, . . . ,k*,k^, . . . , k^_^) 
is a base in Z". Therefore, for the volume of the parallelepiped spanned by 
{k\, . . . ,k*,k^, . . . , k:^_j.), one has the following formula: 

volnikl, ...,k;,ki,..., k^_^) = vol(R"/Z") = 1. 
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On the other hand, since k* — i* € T^"*" for any i = 1, . . . ,r, and the systems 
{£1, . . . ,£*.) and {k^, . . . , A;^_r) are mutually orthogonal, we get 

volniki, . . . ,k^, ki , . . . , k:^_j.) =voln(^x, . . . ,1^., kj^, . . . , k:^_j.) 

=volr(^i, . . . ,^*)vol„_r(A;j'-, . . . , k^_^) 

=voi(y/r*)voi(y^/r^). 

Thus, we have 

^ ' ' voi(y/r*) 

In order to complete the proof of (|3.2p . it remains to apply the well known 
relation 

voi(y/r)voi(y/r*) = i. 

Thus, we can write 

(3.3) n,[S)= ^e{S.l*), 

7*er* 

where 

n,(5,7*) = #(T,(5)nZ!;.). 

Note that, since S is bounded, the sum in the right hand side of (j3.3p has 
finitely many non- vanishing terms. 

Fix 7* G r* . Let xs^* be the indicator of the set S"^* = 5* P| P^* . It is 
easy to see that 

ne[S,-i*) = XT,(S^*){k) = J]] XT, (5^0 (^7* +7) 

= Xs,,{TMkr+l)) 
7eri 

= Y xs.Akr + (re-i(V) - V) + ?;-i(7)))- 

7er^ 

The space V'^ decomposes into the direct sum 

(3.4) v^=Fv^H, 

where Fy = F n V-^. We will write the decomposition of x € V-^, corre- 
sponding to p.4p . as follows: 

X = Xp + Xh, Xf G Fy, Xh € H. 

Note that 

Te(x) = Xp + e~^XH- 

Let p € Cg"(M) be an even function such that < p{x) < 1 for any 
X G M and suppp C (—1, 1). For any tp > and tn > 0, define a function 
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Pt.,tH e C^{V^) by 

where the constant c > is chosen so that Jy± pi,i{x) dx = 1. The function 
Ptp,tH is supported in the eUipsoid 

B{0,tF,tH) = IxeV^ ■.^ + ^ <l\ . 

Define the function ne,tF,tHi^-'l*) by 

where the function XTs{S^,) * PtF,tH ^ ^o^i^'y*) is defined by 

(XTe(5^*) * PtF,tH){y) = J^^ XT,{S^,){y - x)ptF,tH{x) dx, yePj*. 

For any domain S C P^* and for any tp > and tn > 0, denote 

StF,tH = \J{x + BiO,tF,tH)), 
xes 

and 

S-tp-tH = ^7* \ (^7* \ S)tF,tH- 

It is easy to see that, for any e > 0, ti? > and tn > 0, one has 

TeiStF,etH) = {Te{S))tp,tH- 

Lemma 3.1. For any e > 0, tp > and tn > 0, the following inequalities 
hold: 

ne,tp,tH{{Sr)-tF-etH^^*) < rieiS,!*) < ne,tF,tH{{Sj*)tF,etH,l*)- 

Proof. Suppose that A; € n Ts{S^*). For any x G B{0,tF,tH); the point 
k-x belongs to {Te{S^*))tp,tH- Therefore, XTs{{S.^*)tp ,etff)i^ — x) — 1 and 

{XTei{S.,*)tp,,t„) * PtF,tH){k) = 1 Ptp,tH{x)dx = 1 = XT,{S^.){k). 

lik^'LT\T,{S^*), then 

XT,{S^*){k) = < {XTe{{S^*)tp,etH)* PtF^tH){k). 

We get 
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ke 



The second inequality follows of the proven one applied to the domain B\S, 
where S is a sufficiently large ball, containing S. □ 
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For any / € S{V'^), define its Fourier transform / € S{V'^) by 
f{0= f e-2-«'-')/Wdx. 



Recall the Poisson summation formula 

(3.5) Y.f(k) = j^ 

where F"*-* C V'^ is the dual lattice of F"*-, and we used the relation ()3.2p . 
Let us apply (jS.Sp to the function 

(3.6) f{x) = ixT,{(S-,*)tp.etH) * PtF,t„){k^* +x), x£ V^. 

The formula ()3.5p can be applied, because, for any > 0, we have the 
estimate 

(3.7) |^,,,.(0|<C.^^^^^-^^^^, 

One has the relation 

XT.((v (0 = ^-'^e2^^^^'^'-'^^^^'-*^^X(v).„.,(r.(0)- 

Indeed, for any subset S C V-^, we have 

XTAS)iO = e-2-«'-)xT,(5)(V +^)rfx 

'^""'^^'"hsiky* + {T,-i{k-y*) - ky*) + T,-i{x))dx 
-2.^{c,Tax')+(T.(V)-V))x5(^, + x')dx 



e~ 



We also have 



By these relations, the Poisson formula (|3.5p applied to the function / given 
by (|3.6p is written as 



(3.8) ne,tp,tHiiSr)tF,etH,'y*) 

= W E e2-('=*'(^-^^)(V));e(v).,,,,(Te(A:))/)i,i(tpfei. + tH%). 

The series in the right hand side of (j3.8p converges uniformly by the esti- 
mate ([3771) . 
Let us write 



n 



e,tp,tH 



( (^7* )tF ,etH > 7* ) = n',^tF,tH ( (^7* )iF ,etH > 7* ) + <t^,t^ ( (^7* )iF ,etH > 7* ) , 
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where 

= W X{S^*)tp,etuiTeik))pl,l(tFkF+tHkH), 

feer-L*,fe^=o 

and 

Let /c e r-L* be such that % = 0. Then A; G Fy. Since T-^* C Q" and 
n Q" = {0}, we get k = 0. Thus, we have 

'^£,tj.,tjf(('5'T*)tF,£t/f>7*) 

~ |Q| '^('S'-y*)tF,etH 

= -|gj-VOln-r(('S'-y*)tj,,£t^) 

We have the estimate 

Voln-ri{S^*)tF,stH XS-y*) < C{tF + tne), 
therefore, we obtain that 

(3.9) 7*) = |^voU_.(P^* n 5) + OitFS-'i + tHE^-'i). 

Consider the case when k G F-*-* and ku 7^ 0. For any t E Fy and for any 
domain D C P-y* , we denote 

D{t) = {XH e H : k^* + t + XH e D} C H. 

For any function G S{H), denote by Fh{(I)) G its Fourier transform: 

Jh 

It is easy to see that 

[FHiXDit)mH) = I XDikr + t + x,^)e-2'^^(««'^«) dxH, e H. 

J H 
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Therefore, we get 

J Fy J M 
J Fy 

and, as a consequence, we obtain the estimate 

(3.10) \Xis,,)^,,,JUk))\ < [ \FH[xis,,),^,.,„(.,)]{e-'kH)\dxF. 

Hence, our considerations are reduced to the sufficiently well studied prob- 
lem of estimating the Fourier transform of the indicator of a domain, and 
we can apply existing results. 

Proof of Theorem \l.l\ For any sufficiently small e > 0, ti? > and tn > 
0, the domain has smooth boundary. The Stokes formula 

allows us to write the Fourier transform -F/i'[x(5^*)t^ et^C^^F)] oscillating 
integral over the boundary of 



(3.11) FH[Xis.^,),^,,t„{xp)]{-r^) 



1 



r 

that implies the estimate 



e 

d[{S^*)tp,et„{xF)] 



(3.12) |FH[x(w..,{-.)](Oi = o(ier'), lei^^. 

Hence, using the estimates (|3.10p . ()3.12p and ()3.7p . we get 

(3.13) Ktp,t,{{S^')tp.etu,l*)\ 



1 



Putting tn = ch to be a constant, independent of e > 0, and tp = cpe"', 
where 

1 

a = — r, 

p — r + 1 

by p.9p and ()3.13p . we obtain the statement of the theorem. □ 
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Proof of Theorem \1.2[ By assumption, the domain S-y*ixF) = S (1 {7* + 
XF + H} is strictly convex. Therefore, for any sufficiently small e > 0, 
tp > and tn > 0, the domain {S-y*)tF,etHi^F) is strictly convex. By 
the stationary phase method, we derive from the representation (j3.1ip the 
following estimate 



(3.14) = 0(1^1-^'^+')/'), \(\ ^ 00. 

Thus, using the estimate (|3.1Up . (|3.14p and (|3.7p . we obtain that 

(3-15) \nltF,tHiiSr*)tF,etH,7*)\ 

<CF-'i V p{g+i)/2|L„|-{g+i)/2 \ 



N 

dx F dXH 



Put If = e"^, tff = e"^, where Oi? > and an > are chosen as follows. 
If > p — r, then 

2q q-l-2{p-r) 
O-F = — , , , T, Q-H ^ 



g + 1 + 2(p - r) ' q+l + 2{p-r)' 

If 2zJ. < p — r, then 

g + 1 . 
2{p-r + 1) 

Using the estimates (|3.9p and (j3.15p , we immediately conclude the proof. □ 
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